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Abstract. 

We present the evolution of rotation in models of massive single stars covering a wide range of masses and 
metallicities. These models reproduce very well observations during the early stages of the evolution (in particul ar 
WR populations and ratio between type II and type Ib,c at different metallicities, see lMevnet fc Maederll200filh 
Our models predict the production of fast rotating black holes. Models with large initial masses or high metallicity 
end their life with less angular momentum in their central remnant with respect to the break-up limit for the 
remnant. Many WR star models satisfy the three main criteria (black hole formation, loss of hydrogen-rich 
envelope and enough angular momentum t o form an accr etion disk around the black hole) for gamma-ray bursts 
(GRB) production via the collapsar model Considering all types of WR stars as GRB progenitors, 

there would be too many GRBs compare d to observations. If we consider only WO stars (type Ic supernovae as 
is the case for SN2003dh/GRB030329, see lMazzali et alJl2003T) as GRBs progenitors, the GRBs production rates 
are in much better agreement with observations. WO stars are produced only at low metallicities in the present 
grid of models. This prediction can be tested by future observations. 

Key words. : Stars: evolution, rotation, Wolf-Rayet, neutron, black holes: theory Gamma rays: theory, bursts, 
supernova:general 



1. Introduction 

In the previous papers in the series on stellar evolution 
with rotation, the effects of rotation on the stellar evo- 
lution was studied with an emphasis on the early stages 
of the pre-supernova evolution (main sequence, MS, and 
He-burning). The models reproduce very well many ob- 
servational fe atures at various metallic ities, like surface 
enrichments l|Mevnet fc Maederl l2002bl). ratios between 
red and blue supergian ts llMae'derfc^Ievnetll200ll) (here- 
inafter paper VIII . In | M eynet fc Maederl 1|2005|) (paper 
XI) and lMevnet fc Maederl fc003ft (paper XI. the Wolf- 
Rayet (WR hereinafter) star population at different metal- 
licities are studied. In view of the lower mass loss rates ob- 
tained when clumping effects in the winds are accounted 
for, models without rotation do not reproduce the WR 
populations in galaxies. Models of rotating massive stars 
give much better fit to the populations of WR stars at dif- 
ferent metallicities than non-rotating models. They can 
reproduce the observed number ratio of WR to O-type 
stars, the observed ratio of WN to WC stars (for metallic- 
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ities lower than solar), the observed fraction of WR stars 
in the transition phase WN/WC and finally the observed 
ratio of type lb and Ic to type II supernovae. A good fit 
of the observed properties of WR stars is of particular 
interest in this study since WR stars are thought to be 
progenitors of GRBs. 

In lHirschi et alJ l|2004l) (hereinafter paper XII), we de- 
scribed the recent modifications made to the Geneva code 
and the evolution of the models until silicon burning. In 
this paper, we look at the evolution of rotation in massive 
stars with an emphasis on the final stages of the evolu- 
tion towards neutron stars (NS), black holes (BH) and 
especially long soft 7— ray bursts (GRB). 

1.1. Neutrons stars 

There are many observations of neutron stars and pul- 
sars available. A catalogue of observed pulsars an d their 
prope rties is available at the web page of the group IaTNFI 
( 200^) . The fastest young pulsa rs have periods larger than 
10 ms: 16 ms for J0537-6910 llMarshall et al.lll998lh 33 
ms for the Crab pulsar (B0531+21. IStaelin fc Reifensteml 
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Il968l) . 50 ms for J0540-6919 ilSeward et alJll984h . There 
are about 20 pulsars with a period smaller than 100 ms 
in the catalogue with an age less than lOO'OOO years. 
The pulsars may have slightly slowed down since their 
birth but the initial period of pulsars is at least 10 ms. 
NS are supposed to form during the collapse of stars 
with an initial mass between about 10 and 25 Mq. If we 
compare 10 ms with the values obtained in Tables 0H3 
(columns 11 and 12), we see that in general our models 
have much more angular momentum than the observed 
pulsars. The difference can reach a factor of about 100 
around 15 Mq and a factor of about 10 around 60 Mq. 
In order to reconcile the model predictions with the ob- 
servations, additional angular momentum has to be lost 
before the formation of the pulsar. This can occur during 
the pre-supernova stages due to the effects of the mag- 
netic fields not taken into account in this work or during 
the collapse and the explosion. The existence of a primor- 
dial magnetic field in some stars l(Braithwaite fc Spruid 
2004) may also slow down their core. Models including 
the effec ts of magnetic fields have only recently bee n de- 
veloped l|Heger et al.ll2004t iMaeder fe Mevned 12004]) due 
to the complexi ty of the int e rplay between rotation and 
magnetic fields. iHeger et alJ ((2005) show that the brak- 
ing by magnetic fields (which are produced by differential 
rotation during the evolution, see lSpruitll2002l) reduce sig- 
nificantly the discrepancy between the predicted and ob- 
served pulsar periods. Ho wever, these models a s well as 
binary models presented in lPetrovic et alJ l)2005[) contain 
too little angular momentum at the pre-supernova stage 
to produce GRBs. 

In this paper, we study the case without magnetic field, 
which is more favorable to the GRB production. The case 
with magnetic braking will be treated in a further study. 
In addition, the present study offers a basis of comparison 
for future models including magnetic fields. It is impor- 
tant here to note that internal magnetic braking is most 
efficient during the Red SuperGiant (RSG) phase during 
which the star has a fast rotating core and a slow rotat- 
ing envelope. This means that magnetic braking is much 
less effect ive for stars which s kip the RSG stage, as men- 
tioned in IHeger et al.l l|2005() . The stars which skip the 
RSG stage are the most massive WR stars, with masses 
larger than about 60 M@. Furthermore, more massive 
stars have shorter lifetimes and thus less time to be slowed 
down. It is also possible that mechanisms that slow down 
the nascent neutron star may not be as efficient when a 
black hole is formed. These differences between stars with 
masses around 15 Mq (forming a neutron star and going 
through the RSG stage) and 60 Mq (not going through 
the RSG phase and forming black holes) could possibly ex- 
plain, with the same physics, the pulsar periods and the 
GRBs. Further studies will verify this possibility. 

Braking between the collapse and the pulsar formation 
can also occur via different mechanisms: r-Mode insta- 
bility, neutrino-powered magnetic stellar winds, fall-back 
and the propeller mechanism. However, these mechanisms 
may not be efficient enough to slow down the core ef- 



ficiently (see lHeger et al.l l2005. for a discussion). Latest 
models studying gravitational waves seem to indi cate that 
braki ng is possible during and after the collapse l)Ott et alJ 
2005). If the cores are not slowed down during the pre- 
supernova stages, rotation may play an important role in 
supernova explosions. Rotation may in this case provide a 
large amount of energy for the explosion [~ 10 52 erg) but 



this is generally not observed Uanka et alJl2005j) . 



1.2. Black holes and GRBs 

Theoretical m odels still struggle to reproduce supernova 
explosions fsee lFrver fe Hegerll2005l for a discussion and 
references). It is therefore not possible for us to predict 
with certainty or precision the fate of the remnant of our 
mode l s. Nev erthe less, following the e stimates presented in 
iFrverl lll999l) and lHeger etllfll2003h . we will consider in 
this study that the lower mass limit for black hole for- 
mation is around 25 Mq. The upper mass limit depends 
strongly on the mass loss and is around 100 Mq at solar 
metallicity. T he maximum neutron star mass is also still 
uncertain (see |Sriniyasa n1 l2002t iMorrison et"aTH 2004. for 
recent studies). Depending on the nuclear equation of state 
used and rotation rate, the upper mass limit for neutron 
stars can vary between 1.6 and 3 Mq although 2 to 2.5 Mq 
is a more common range. Black holes cannot be observed 
directly. Thus they are usually detected in X-ray binaries. 
In certain cases, observations allow the determination of 
the masses of the binary companions. If the accreting ob- 
ject has a derived mass larger than 2-3 Mq, it is consid- 
ered as a BH candidate. T he best known candidates are 
LMC X-3 and Cyg X-l Csee lKaoer fc van der Meerll2005L 
and references therein). 

GRBs can be divided in two m ain types: a) short 
and hard and b) long and soft. See IPiranl l(2005h for a 
recent review on GRBs. The long soft GRBs have re- 
cently been connected to super nova (SN) explo sions of 
the type Ib,c (see for example iMathesonl |2003|) . Since 
then, several studies h aye been devoted t o find which 
stars can le ad to GRBs. IHeger et all l|2004^ . IHeger et alJ 
( 200sJ) and iHirschi et al.1 l|2004f) looked at massive sin- 
gle star at solar metallicity as progenitors. Note that the 
two models differ in the treatment o f meridional circu- 
lation, which is an advective process ijMevnet fe Maederl 
2002a). Our models account for the advection of angu- 
lar momentum during the MS phase while models of 
Heger and collabor ators treat meridional cir culation as 
a diffusive process. IMaeder fe Mevnetl l)2005[) show that 
the treatment of meridional circulation as an advective 
process is crucial to model the interplay between circu- 
lation and magnetic braking. Another difference in this 
article is the fact that we study in detail the effects of 
metallicity in the context of GRBs for models of differ- 
entially rotating stars. This is important since metallici- 
ties lower than solar are expected to be more favourable 
to the production of GRBs (due to weaker mass loss , 
seelMa,cFadven fc WooslevirSlWooslev fc He g erll2003h . 
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IPodsiadlowski et all l|2004|) and llzzard et all l)2004|) con- 
sidered both single and binary massive stars. Note that 
llzzard et alJ l|2004j) assume that the loss of angular mo- 
mentum is proportional to the amount of mass lost and 
do not consider internal transport of angular momentum. 
Their approximate treatment of the evolution of angular 
momentum leads them to the conclusion that only binary 
stars can retain enough angular momentum in their core to 
form GRBs. Even without magnetic braking, their mod- 
els predict no GRBs from single massive star s which is in 
contra d iction with our mode ls and those of iHeeer et alJ 
l|2000|) . IPetrovic et"aTl l)2005|) look at single and binary 
systems with and without magnetic fields at solar metal- 
licity. They find that in the models followed, both single 
and binary models without magnetic field can produce 
GRBs and both single and binary models with magnetic 
field considered in the study cannot produce GRBs. They 
conclude that, if the present modelling of magnetic fields is 
accurate, GRBs have to be produced in s ome e xotic chan- 
nels of binary systems. iFrver Hegel J20O5) looked at 
massive binary star mergers, which are believed to be one 
of the best binary candidates for GRBs. They follow the 
evolu tion prior to and a fter the merger with the KEPLER 
code llHeeer et al-l fcoOO) and the merger process with 3D 
simulations, which is very interesting and which has a lot 
of potential for further applications. They find that, in 
some cases, merged helium stars can retain more angular 
momentum than single massive stars. However, they ob- 
tain these faster rotators only for models in which they 
remove angular momentum artificially without removing 
mass from the merged star. If more reasonable mass loss 
prescriptions are used after the merger, the final angular 
momentum contained in the central remnant is similar in 
merged stars and in single massive stars. This is due to 
the strong mass loss that WR stars experience, especially 
the fast rotating ones. Merged stars lose a large amount 
of mass and some merged systems probably form NSs in- 
stead of BHs. It is therefore not proven that binary stars 
retain more angular momentum than single stars and the 
exact frequency of potential binary star progenitors has 
not been compared quantitatively with the GRB rates. 

Here we discuss the possibility for single massive stars 
to be progenitors studying models covering a large range 
of initial masses and metallicities. These models repro- 
duce very well the WR star population and many of them 
retain enough angular momentum to produce GRBs as 
shown below. We described the calculations in Sect. 2, the 
evolution of rotation in Sect. 3 and present our predicted 
GRB rates in Sect. 4. 

2. Description of the calculations 

In order to study the evolution of core rotation and its 
dependence on initi al mass and metall i city, w e used the 
mo dels presented in IMevnet fc Maederl (|2003j) (paper X) 
and lHirschi et all ll2004l ) (pap er XII) for solar metallicity 
and iMevnet fc Maedeiri|2005h (paper XI) for non-solar 
metallicities, which span a wide range in initial masses 



(9 to 120 Mq at solar metallicity) and metallicities (Z = 
0.004 to Z = 0.040). 

All these models (except for one 60 Mq model 
at Z = 0.004) have an initial rotational velocity of 
300kms _1 , which gives an average value on the MS of 
about 220 km s -1 , the average observed rota tional veloci- 
ties on the MS (see for example iFukudalU 982h . They were 
all calculated with the same input physics: Schwarzschild 
criterion for convection, overshooting of 0.1 Hp for the hy- 
drogen (H) and helium (He) burning cores, same mass loss 
prescriptions, treatment of rotation including meridional 
circulation and shear diffusion (see papers X-XII for more 
details) . 

Non-solar metallicity models have been evolved until 
the end of core He-burning. Solar metallicity models have 
been evolved until the end of core He-burning for the 9, 
85 and 120 Mq models, end of oxygen (O) burning for 
the 12 Mq model and end of silicon (Si) burning for the 
15, 20, 25, 40 and 60 Mq models. As we shall see later, 
the largest changes in the angular momentum of the core 
occurs during H and He-burnings and this is why we did 
not evolve non-solar metallicity models further than the 
end of He-burning. 

Several quantities are useful for studying the evolution 
of rotation: 

- The angular velocity at mass coordinate m or ra- 
dius r (here r is the average radius of a shell, see 
IMevnet fc Maederl 19"97i for details), £1 [s -1 ], and its ratio 
to the Keplerian velocity, $1/SIk, where VI k = \/Gmjr i . 
When Q,/VLk becomes larger than about 0.9, the star gets 
elongated along its equator and at break-up, the star 
equatorial radius is 1.5 times larger than the polar ra- 
dius (value obtained using the Roche model). Therefore 
the star may reach break-up before SI/Qk approaches 
one, in particular for models with a large Eddington fac- 
tor which reach break-u p for much lower than fix 
l|Maeder fc MevnedExxifr . 

- The angular momentum of a core of mass M, Cm = 
f jmdm [gcm 2 s _1 ], where j m = 2/3fi m r^ [cm 2 s^ 1 ] 
is the specific angular momentum at mass coordinate m. 
The momentum of inertia of a core of mass M, Tm = 
J^2/3r 2 dm [gem 2 ]. 

The properties of the models are presented in four 
Tables: Table0]for the metallicity of the Small Magellanic 
Cloud (SMC, Z=0.004), TableElfor the Large Magellanic 
Cloud (LMC, Z=0.008), Table for solar metallicity 
(Z=0.020) and Table □ for the Galactic center (GC, 
Z=0.040). In these tables, for each model, we give the 
initial mass and velocity as well as the remnant mass es- 
timated from th e carbon-oxyg en (CO) core mass, using 
the relation from lMaederl {.1992). We used the value of the 
CO core at the end of the evolution for the calculation (at 
the end of core He-burning in general and at the end of 
Si-burning for the models from paper XII). Then, 

— column 1 is the evolutionary stage to which the values 
correspond, 

— column 2 is the total mass at the given stage, 
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Remnant mass (M co ), D ini = 300 km s 1 
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Fig. 1. Remnant gravitational masses (estimated from 
the CO core masses) as a function of the initial mass of 
the models. 

— column 3 is the angular momentum contained in 
the remnant, £ rcm = J^ Ilcln 2 / 3 fl r 2 dm in units of 
10 50 [gcn^s" 1 ], 

— column 4 is the moment of inertia of the remnant, 
I rcm = J M " m 2/3 r 2 dm in units of 10 55 [gem 2 ], 

— column 5 is the specific angular momentum at the rem- 
nant edge, jrom, in units of 10 16 [cm 2 s _1 ], 

— column 6 is the average specific angular momen- 
tum in the remnant, j rom = C lem /M Tem , in units of 
10 16 [cir^s -1 ], 

— column 7 is the specific angular momentum at the mass 
coordinate 1.56 M©, ji.56, in units of 10 16 [cm s ], 

— column 8 is the angular momentum contained in the 
inner 1.56 M , £i. 56 = f^ me 2/3 fir 2 dm in units 
of 10 50 [gcm 2 s- 1 ], 

— column 9 is the specific angular momentum at the mass 
coordinate 2.5 M@, j'2.5, in units of 10 16 [cm s ], 

— column 10 is the angular momentum contained in the 
inner 2.5 M Q , £2.5 = Jq M<s 2/3 fir 2 dm in units of 
10 50 [gcm 2 s- 1 ]. 

Assuming that a neutron star with a baryonic mass, M& = 
1.56Mq, and with a radius, R = 12km, would form from 
the models, we calculated the two following quantities: 

— (column 11) the ratio of the NS angular velocity to its 
Keplerian angular velocity, O/Qjf (NS). 

— (column 12) the initial rotation period of the neutron 
star, V ro t in units of milli seconds. 

A proto-neutron star with a baryonic mass, Mb = 
1.56 Mq, loses a binding energy (BE) of 0.159 M Q due to 
neutri nos emission. This va l ue is c alculated using Eq. (36) 
from ILattimer fc Prakashl (|200ltt : BE/M = 0.6/3/(1 - 



Table 1. Baryonic masses, Mb, gravitational masses, M g 
and the BE (all in M Q units, M Q c 2 = 1.79 1 54 erg) of 
the different models calculated using Eq. Q and R = 
15.12 M b _1 ^ 3 (Newtonian Dolvtropes lShapiro fc Teukolskvl 
11 9831 p. 245). 



A/ini M b M g BE 



Z=0.004 



30 


3.7284 


2.8140 


0.9144 


40 


5.3954 


3.5430 


1.8524 


60 


8 5150 


4 3964 


4 1186 


60" 


3.7916 


2.8464 


0.9452 


120 


5.1429 


3.4478 


1.6951 




Z= 


0.008 




30 


3.7256 


2.8125 


0.9130 


40 


5.2110 


3.4740 


1.7370 


60 


4.9085 


3.3550 


1.5535 


120 


4.0571 


2.9780 


1.0792 




Z= 


0.020 




9 


1.1808 


1.1012 


0796 


12 


1.4616 


1.3338 


1278 


IS 


1.8496 


1.6361 


0.2135 


20 


2.5661 


2.1382 


0.4280 


25 


3.0558 


2.4424 


0.6135 


40 


3.8527 


2.8773 


0.9754 


60 


4.3228 


3.1027 


1.2201 


85 


3.7764 


2.8386 


0.9377 


120 


3.5392 


2.7145 


0.8247 




Z= 


0.040 




20 


2.6804 


2.2119 


0.4685 


25 


3.1257 


2.4834 


0.6423 


40 


3.5718 


2.7320 


0.8399 


85 


2.5705 


2.1410 


0.4295 


120 


2.5329 


2.1164 


0.4165 




Z=0.040, Mwr(2) 




40 


3.0225 


2.4226 


0.5999 


60 


1.9391 


1.7027 


0.2364 


85 


1.9520 


1.7122 


0.2398 


120 


1.9637 


1.7208 


0.2429 




with Ui n 


= 500 km s 


-1 



0.5/3), in which M = M g corresponds to the gravitational 
mass and where (3 = G M g /Rc 2 , and using the fact that: 
Mb = BE + M g . These equations give one second degree 
equation: 

0.1 G M 2 -{Re 2 + 0.5 G M b )M g -Rc 2 M b = (1) 

from which we can find the gravitational mass, M g = 
I.4OIM0, which is very close to the average obser- 
yational value of neutron s tar (gravitational) masses 
(Kap er fc van derMee][2005). This is why we chose the 
value of 1.56 for the baryonic mas s. The radius of 12 
km wa s estimated using Fig. 2 from ILattimer fc Prakashl 
l|200lh . The NS momentum of inertia is calculated us- 
ing Eq. (29) from ILattimer fc PrakasbJ l|200ll) . which 
for the mass and radius chosen here corresponds to 
Jns(TVII)/M R 2 ~ 0.36. Since the core loses mass dur- 
ing the collapse, it also loses angular momentum. We 
assumed here that the angular momentum loss due to 
neutrino losses is proportional to the quantity of mass 
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Fig. 2. Evolution of the angular velocity as a function of the radius at different stages of the evolution. Left: the 15 
M Q model becomes and ends up as a RSG. It therefore has a strong differential rotation between the core and the 
envelope. Right: the 120 Mq model becomes a WR before the end of the MS and will not go through the RSG stage. 
This last model does not have a strong differential rotation between the core and the envelope and should not be 
slowed down significantly by magnetic fields. 



lost. Therefore the angular momentum of the NS, £ns — 
£i.56 {Mg/Mb). This assumes that neutrinos are emitted 
uniformly from the entire core. This is probably correct 
but additional angular momentum may be lost due to the 
inter action of these neu trinos with the outer layers of the 
core llHee-er et a,l.l|2005|) . 

Assuming uniform rotation in the neutron star and us- 
ing the relation £ N s = Ins ^ns > we have fl^s = £ns / ^ NS ■ 
Taking into account no angular momentum losses other 
than those due to neutrinos described above, 

Ons=A.56 (M g /M b )/(0M R 2 M g )=A. 56 /(0-36 R 2 M b ). 



The Keplerian angular velocity is 17k = \J G M g /R 3 . 
Finally, the initial period of rotation of the neutron star 
is given by V ro t = 2 7t/^ns- 

We also give values of the angular momentum at and in 
the inner 2.5 Mq (near maximum mass for neutron stars). 
These results allow to compare different models at a given 
mass coordinate rather than at the remnant edges which 
span a wide range of masses (see Tables 0H?|). Finally, we 
used Eq. to obtain the gravitational masses, M g and 
the BE (all in Mq units, Mq c 2 = 1.79 10 54 erg) for our 

— 1/3 

models using R = 15.12 M u (Newtonian polytropes 
IShapiro fc Teukolskvll98l p. 245) These values are given 
in Table^ an d the gravitational masses are shown in Fig. 
n We see that, using our estimates, black holes are ex- 
pected to form from stars with masses above about 20 Mq 
for solar or lower metallicities. At higher metallicities, the 



formation of a black hole or neutron star depends on the 
mass loss prescription. At Z=0.040 and using a mass loss 
dependent on metallicity (models MZ), we expect black 
hole to form from stars with masses between about 20 
and 55 M e . 

3. Evolution of rotation 

In paper XII, we presented the evolution of the angular 
velocity profiles through the pre-supernova stages for the 
solar metallicity 25 Mq model. We showed that the an- 
gular velocity increases in the core via successive contrac- 
tions. In the outer part of the star, after the MS, f2 either 
has very small values, when the star is a red supergiant 
(RSG), or has values similar to the initial values when 
the H-rich envelope of the star is lost and the star be- 
comes a WR star. If a star is very massive, the star may 
become a WR before the end of the main sequence (MS) 
and skip the RSG phase (see Fig.FJJ. This point has im- 

of st ars by magnetic 
20051) . Indeed in the 
theory elaborated bv ISpruitl 1)2002|) . differential rotation 
produces strong magnetic fields, which in turn slow down 
the core of the star. This means that stars that go through 
the RSG phase will experience a strong braking due to the 
large differential rotation between the core and the enve- 
lope while stars that skip the RSG stage might be slowed 
down much less. 

The evolution of the angular momentum profiles has 
also been discussed in paper XII. In our models, we con- 
sider three transport processes: convection, shear diffusion 



portant consequences f or the brakin . 
fields, as mentioned in Heger et al.l 



6 



R. Hirschi et al.: Stellar evolution with rotation XIII: 




Fig. 3. Profiles of the specific angular momentum, j r , as a function of the Lagrangian mass coordinate, m, for different 
models: 15 Af© (solid lines), 25 M© (dashed lines) and 60 M (dotted lines). Left: On the ZAMS. Right: At the end of 
Si-burning, js — \fY2.Gmj c (dashed-dotted line) is the specific angular momentum necessary for matter to form an 
accretion disk around a non-rotating black hole (Schwarzschild metric) and is displayed here as a reference profile. The 
specific angular momentum varies significantly with the mass coordinate at the end of Si-burning. For the average 
profile, the heavier the initial mass of the star, the stronger the mass loss and angular momentum loss. Therefore 
the core of the 60 M 9 model ends with much less angular momentum than on the ZAMS and loses more angular 
momentum in comparison with lighter models (15, 25 Mq). 



and meridional circulation. Convection makes Q constant 
and therefore transports angular momentum from the in- 
ner part of a convective zone to the outer part of the 
same convective zone. Shear diffusion reduces differential 
rotation and also transports angular momentum outwards. 
Meridional circulation can transport angular momentum 
inwards or outwards. Mass loss acts somewhat indirectly 
on the angular momentum of the core by affecting f2 and 
the gradients of inside the star and thus the efficiency of 
the transport processes. The main conclusions about the 
evolution of the angular momentum are the following: 

- The angular momentum generally decreases through- 
out the pre-supernova evolution. 

- The largest decrease occurs during the MS. 

- After He-burning, only convective zones reshape the 
angular momentum profile and produce teeth along its 
profile (redistributing angular momentum from the inner 
boundary of a convective zone to the outer boundary of 
the same zone) without removing significant amount of 
angular momentum from the core. This shows that the 
angular momentum at the end of He-burning is a good 
approximation of the angular momentum contained in the 
core at the collapse. 



3.1. Dependence on the initial mass 

The variation of both the remnant mass and the specific 
angular momentum as a function of the mass coordinate 
makes the comparison between different models quite dif- 
ficult at the pre-supernova stage (see Fig.[3J). Indeed, the 
specific angular momentum, j, varies significantly in the 
core (see Fig. 7 in paper XII). Until the end of He-burning, 
j increases monotonically with the mass coordinate. At 
the end of Si-burning, the teeth produced by the convec- 
tive zones make the final j at a given mass coordinate os- 
cillate between values slightly above or a few times under 
the values at the end of He-burning. The remnant mass 
also obviously varies with the initial mass (see values in 
Tables HQ. 

This is why we also give the values of the specific an- 
gular momentum at fixed Lagrangian mass coordinates 
(1.56 and 2.5) and give the average value inside the rem- 
nant (j rcm ). It is easier to compare between models at the 
end of He-burning and this is what we do in the next 
section. Nevertheless, we present the values at the end of 
Si-burning for solar mctallicity models to show how the 
specific angular momentum at a given mass coordinate 
can vary during the advanced stages. 

The general dependence on the initial mass of the final 
specific angular momentum in the core can be seen in Fig. 
[3J The larger the initial mass, the smaller the final specific 
angular momentum in the core. This is due to a stronger 



R. Hirschi et al.: Stellar evolution with rotation XIII: 



7 



M= 40 M m , u. , = 300 kms" 1 , a„„ =0.1 




Z = 0.004 ZAMS 

Z = 0.040 ZAMS 

Z = 0.004 end He 

Z = 0.040 end He 



10 15 

m r [MJ 

Fig. 4. Profiles of the specific angular momentum, j r , as 
a function of the Lagrangian mass coordinate, m, at the 
ZAMS and at the end of He-burning for 40 M models at 
Z=0.004 and Z=0.040. The metallicity dependence of the 
mass loss affects the evolution of the profile, especially at 
the surface. The decrease of angular momentum is larger 
at higher metallicity due to a larger mass loss and a more 
efficient transport. The specific angular momentum neces- 
sary to form a accretion disk around a non-rotating black 
hole (Schwarzschild metric), jg = \/V2GM/c is given for 
reference. We see that a low metallicity star is more likely 
to form an accretion disk than a high metallicity star. 

mass loss and a more efficient transport of angular momen- 
tum out of the core. The more efficient transport in more 
massive stars can be explained by the smaller compactness 
of these stars. Indeed, higher compactness reduces the out- 
ward transport by meridional circulation (Gratton-Opik 
term proportional to see paper VII). 

3.2. Dependence on the metallicity 

Figure 0] shows the initial and final angular momentum 
profiles (as a function of the Lagrangian mass coordinate) 
for 40 M Q models at metallicities equal to 0.004 (SMC) 
and 0.040 (Galactic center). Since higher metallicity im- 
plies stronger mass loss and less compact stars, the depen- 
dence of the final specific angular momentum in the core 
on the metallicity is similar to the dependence on the ini- 
tial mass of the star. Therefore, the higher the metallicity, 
the smaller the final angular momentum in the core. 

4. Models at the final stages and predictions for 
long soft gamma ray bursts (GRBs) 

Stellar mass black hole candidates (Cyg X-3, LMC 
X-l, ...) are mostly observed in X-ray binaries (see 



iKaper fc van der Meerl 120051 and references therein). 
Usually, a lower limit for the mass of the accreting body 
which is larger than about 2 to 3 solar masses makes it 
a black hole candidate. Observationally derived masses 
range from about 3 to 15 (for GRS1915+105). The gravi- 
tational masses in Fig. ^ (and Table can give us a very 
rough idea of which models may form BHs. Black holes 
are expected to form from stars with masses above about 
20 Mq. The upper mass limit for BH formation depends 
on the mass loss prescription. At Z=0.040 and using a 
mass loss dependent on metallicity (models MZ), we ex- 
pect black hole to form from stars with masses between 
about 20 and 55 Mq (consdering models with remnant 
masses larger than 2 M in Fig.rj). 

No observation of BH rotation is available. 
Theoretically, one can use the Kerr metric to de- 
scribe rotating black holes. In this metric, the parameter 
a = Jc/GM denotes the angular momentum per 
unit mass (in units of gram if the other quantities are 
in CGS units: G, the gravitational constant, equals 
6.67 10~ 8 cm 3 s -2 g -1 , c, the speed of light, equals 
2.9979 10 10 cms~ 1 , M, the remnant mass, is in grams 
and J, the angular momentum of the remnant is in 
gcm 2 s _1 ). In the Kerr metric, a/M can take values 
between (non-rotating BH: the Kerr metric becomes 
equivalent to the Schwarzschild metric) and 1 (1 for 
maximally rotating BHs). Figure [S] (right) shows the 
predicted values for a/M, assuming no loss of angular 
momentum during collapse (apart from a uniform loss 
due to neutrino emission). Stars lighter than about 20 Mq 
probably form neutron stars. We see that models between 
20 and 60 M have values exceeding the maximal value. 
These models probably form an accretion disk before the 
black hole is fully formed. 

Accretion disks around BHs are very interesting 
in the context of long soft GRBs. These bursts 
were recently c onnected with SNe of type Ib,c (see 
iMatheso n] |2003l for exam ple). Three current scenar- 
ios dWooslev fc Hegerll2003h for GRB production are the 
supranovae model llVietri fc Stella! 1 1999J) . the magnetar 
model JWheeler et aL 120021) and the collapsar mechanism 
l)Wooslevlll993j) . In this last mechanism, a star collapses 
into a black hole and an accretion disk due to the high 
angular momentum of the core. Accretion from the disk 
onto the central black hole produces bi-polar jets. These 
jets can only reach the surface of the star (and be de- 
tected) if the star loses its hydrogen rich envelope before 
the collapse. WR stars are therefore good candidates for 
collapsar progenitors since they lose their hydrogen rich 
envelope during the pre-SN evolution. Which WR stars 
form BHs and which have enough angular momentum to 
form an accretion disk around them? The first part of 
the question can be answered by looking at the estimated 
gravitational mass of the remnants presented in Fig.^and 
given m Tabled In this study, using our rough estimate for 
the remnant mass, the minimum initial mass for black hole 
formation is between 20 M Q (2 Mq as the maximum mass 
for neutron stars) and around 30 M (2.5 Mq). We use 
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Fig. 5. Left: Specific angular momentum at the edge of the remnant as a function of the initial mass of the models. 
The values are determined from the models at the end of He-burning. For solar metallicity models, values determined 
at the end of Si-burning (last) are also given. Short-long dashed lines represent models w ith Z=0.040 with a mass loss 
during the WR phase which depends on the metallicity (MZ, see lMevnet fc Maederl200.'^ . At the edge of the remnant, 
values at the end of Si-burning are usually lower than at the end of He burning but they can also be slightly larger 
as is the case for the 20 Mq model here. Right: Angular momentum per unit mass, a = J c/GM, divided by the mass 
of the remnant, M, as a function of the initial mass of the models. The horizontal line with a/M — 1 corresponds to 
the maximally rotating BHs. 



for GRB rate predictions the lower limit of 20 Mq. Note 
for the following discussion that the present WO models 
definitely form BHs. In order to form an accretion disk 
around a BH, the matter on the last stable orbit (LSO) 
must have an angular mome ntum larger than jx — ?~lso c 
l|Shapiro fc Teukolskvl Il98i p. 428). The radius of the 
last s table orbit, n.sn, is gi v en by r ms in formula (12.7.24) 
from lShaniro fc Teukolskvl l| 19831 p. 362) for circular or- 
bit in the Kerr metric, tlso depends on a and the direc- 
tion of the orbit. We consider here direct (co-rotating) 
orbits. In this case tlso varies from GM/c 2 for a = 1 
to 5 GM/c 2 for a — 0. Note that, for the Schwarzschild 
metric, r LSO = \/l2 GM/c 2 and j s = VUGM/c. For 
models with a > M, we used tlso = GM/c 2 . Figure 
(left) shows the specific angular momentum at the edge 
of the remnant, j rcm , for the different models. The values 
are derived from models at the end of He-burning. For 
solar metallicity models, values determined at the end of 
Si-burning are also given. At the edge of the remnant, the 
specific angular momentum at the end of Si-burning is 
usually smaller than at the end of He-burning but it can 
also be slightly larger as is the case for the 20 M© model 
here. We therefore use values at the end of He-burning for 
our discussion here. 

Between 15 and about 40 Mq, j rem increases (see Fig. 
|5j). This is due to the increase in the remnant size and 
to the fact that the specific angular momentum increases 
with mass coordinate (see Fig. UJ. If we look at the same 



mass range in Fig. [f)| we see that the lower mass stars ro- 
tate faster with respect to break-up rotation than heavier 
stars. In fact, the trend is monotonic: the heavier the star, 
the slower (compared to the maximally rotating BHs) it 
rotates. This is true not only at the remnant edge but 
also at different mass coordinates: 1.56 Mq (see Tables 
SHU and 2.5 Mq (see Fig. bottom-left). Even though 
this trend is already present on the ZAMS (see Fig.Eltop- 
left), it is much weaker and the evolution is responsible for 
the final results. As said earlier, the stronger mass loss and 
more efficient transport of angular momentum are respon- 
sible for the lower final angular momentum in the core at 
high mctallicities. For the metallicity dependence, Fig. H3 
(top-left) shows that lower metallicity models with a same 
initial mass start with a smaller ratio j/jk (and end with a 
higher one). The lower initial ratio is mainly due to the dif- 
ferent (higher) remnant mass. Indeed, j'k is proportional 
to the remnant mass and j = fir 2 is roughly proportional 
to Mrem (if and p are constant), which means that the 

— 1/3 

ratio j/jx. is roughly proportional to M rC m ■ Lower metal- 
licity models having larger remnant masses start with a 
smaller ratio. In this study we choose the same initial sur- 
face velocity for models of different compactness (lower 
metallicity = higher compactness) and not the same an- 
gular momentum or same ratios j/jk or Q/tt c . Future 
studies will consider models with similar angular momen- 
tum or ratios as well as identical velocity, which is very 
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Fig. 6. Top left: Initial ratio of the specific angular momentum at the edge of the remnant to jx = ^lso c (the value 
of angular momentum necessary to form an accretion disk at the edge of the remnant) on the ZAMS as a function of 
the initial mass of the models. The initial ratio is shown to see whether trends observed at the end are due to initial 
conditions or differences in the evolution. Top right: Final ratio of the specific angular momentum at the edge of the 
remnant to jjc = tlso c as a function of the initial mass of the models. The values are determined from models at 
the end of He-burning. At solar metallicity, values at the end of Si-burning (last) are also given. Bottom left: Ratio 
of the specific angular momentum at the Lagrangian mass coordinate of 2.5 Mq to jx — ^lso c as a function of the 
initial mass of the models. This graph is used to compare models at a fixed Lagrangian mass coordinate. Bottom right: 
Ratio of the specific angular momentum at the edge of the remnant to js — yl2 G M /c as a function of the initial 
mass of the models. In this j/js plot, the maximally rotating BH value of jx/js = 1 / \/T2 is also given for reference. 
Short-long dashed l ines represent models wit h Z=0.040 with a mass loss during the WR phase which depends on the 
metallicity (MZ, see lMevnet & Maederl l2005). At the edge of the remnant, values at the end of Si-burning are usually 
lower than at the end of He-burning but they can also be slightly larger as is the case for the 20 Mq model here. 



important for extremely low metallicity stars. Even if high a smaller one. Therefore, the trend of higher metallicity 
metallicity models begin with a higher ratio, they end with 
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Table 2. Predicted rates [yr -1 ] of WR stars and of long 
soft gamma ray bursts (GRBs) from single star progeni- 
tors as well as limiting masses at various metallicities. 



gle stars. For reference, we use the following SN rate: 





Reference rate: TZ^^^ — 


710~ 3 






ZsMC 


Zlmc 


Zq 


Zgc, 




1.31E-03 
6.33E-04 


1.90E-03 
7.58E-04 


2.29E-03 



2.45E-03 



MSb(WR) 
Kqrb(WR) 


32 
95 

1.15E-03 


25 
95 

1.74E-03 


22 
75 

2.01E-03 


21 

55 

1.92E-03 


M max (WO) 


50 
110 
6.33E-04 


45 
100 
7.58E-04 








Kghb(WO) 


50 
95 

4.74E-04 


45 
95 

5.99E-04 








ijOBS 
'^GRB — 


310" 6 - 610 


-4 a 





a Depending on the beaming angle of GRBs. 

Table 3. Final WR type for WC/WO stars at various 
metallicities. At Zqc, models with metallicity dependent 
mass loss rates were used. 





ZSMC 


2lmc 


Zq 


Zgc 


25/30 M 




WC6 




WC5 


40 Mq 




WC4 


WC4 


WC4 


60 Mq 


WO 


WO 


WC4 


WC7 


85 Mq 




? b 


WC4 


WC6 


120 Mq 


WC4 


WC4 


WC4 


WC7 



means the model docs not become a WC or WO star, 
means the model was not calculated. 



leading to smaller final angular momentum (compared to 
the break-up value) is confirmed. 

Concerning the criterion of angular momentum for 
GRB production, our models show that, the lower the 
initial mass, the more probable the GRB production. The 
specific angular momentum necessary to form an accre- 
tion disk around a BH depends on the rotation of the BH 
itself. If we want to remove this dependence, we can com- 
pare jrem with the Schwarzschild metric js = vl2 GM/c 
(see Fig. bottom-right) although a rotating collapsing 
core must form a rotating black hole. In this j /js plot, the 
maximally rotating BH value of jx/js = l/V 7 !^ is given 
for reference. This plot shows that models below about 40 
Mq even have j/js values larger than one. 

4.1. Predicted rates of long soft gamma ray bursts 
(GRBs) from single star progenitors 

Following the method used in lPodsiadlowski et al.1 l|2004|) . 
we present, in Table |21 predictions for the rate of long 
soft GRBs produced in the collapse of massive sin- 



TpOBS 



CSFR /j 



15OM 
10M o 
-1.351 



-2.35 



dm 



710" 3 yr _1 (= 

CSFR ri50" l a;j - / _ The lue f ? 1Q -3 yr -l 

is the same as in ijPodsiadlowski et"al"1l2004() and repre- 
sents the total number of core collapse SNe per year for an 
average galaxy (the Milky Way galaxy has a slightly higher 
value ~ 1.2 10~ 2 yr -1 , see lCappellaro et ai1 ll999i. csfr is 
a normalisation factor taking into account the star forma- 
tion rate (SFR). It is d etermined by the equation above. 
We use Salpeter's IMF (|Salpeterll955|) and an up per mass 
limit of 150 M© following the recent results (see iKroupal 
200$l and references therein). The other rates are then cal- 
ibrated to this reference rate: for example, 7^grb(WR)= 



(710- 3 //, 



15OM 
10M o 



m 



-2.35 



dm) ■ J } 



MS 



,(WR) 
,(WR) 



' 2 - 35 dm. The 



rates of WR stars are deri ved using the minimum m ass 
limits given in Table 2 from lMevnet fc Maed"erl l)2005|) . 

These minimum masses correspond approximately to 
the minimum mass for GRB production since these stars 
have lost their hydrogen envelope by definition and proba- 
bly form BHs (have remnant gravitational masses > 2 Mq; 
see Fig. ^ and Table 0) . T^grb ( WR) corresponds to the 
rate of GRBs produced by massive single stars which form 
BHs, lose their hydrogen rich envelope and have enough 
angular momentum at the edge of the remnant to form an 
accretion disk. We can see that the rates are much higher 
(10-1000 times) than the observed long soft GRB rate, 
ftgg| ~ 3 10~ 6 - 6 10" 4 llPodsiadlowski et allEool . for 
beaming angles between 15 and 1°. 

An additional constraint for GRB production might 
be that the star must be a WR of type WO at the time 
of the explosion. WO stars are WR stars which have lost 
not only their hydrogen rich envelope but also the he- 
lium rich layers. They explode as type Ic supernovae. 
Although there is no firm evidence to support this fact, 
the link between GRBs and supernovae wa s firmly estab- 
lished for the type Ic hypernova SN2003dh l|Mazzali et alJ 
2003). Furthermore, WO stars have lost both their hy- 
drogen and helium rich layers. This means that their en- 
velo pe is more tenuous, thu s favoring the escape of the 
jets. ISmith k, Maederl l|l99lf) studied the metallicity de- 
pendence of WR stars of type WC and WO. They found 
that the lower the initial metallicity Z, the higher the 
(C+0)/He number ratio when the star first becomes a 
WC star. This is due to the fact that at low metallicity, 
mass loss is smaller and the newly synthesised C and O are 
revealed later in the evolution. It implies that WO stars 
are mainly produced at low metallicity. In the present 
grid of models, no WO star appears at solar or higher 
metallicities (see Tableland Fig. 01. In Fig. H we can 
roughly estimate the mass ranges of WO star formation 
at the different metallicities. The rate of WO stars and the 
mass ranges are given in Table |5| Combining all the cri- 
teria for GRB progenitors, we obtain the predicted rates 
for GRB production in a normal galaxy coming from WO 
stars, ^grb(WO) given in Table El 

These rates are upper limits because we did not con- 
volve our rates with the initial rotational velocity distribu- 
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Fig. 8. Final surface(C+0)/He number ratio as a func- 
tion of the initial mass of the models at different metallic- 
ities. WO stars are only found at the metallicities of the 
Magellanic clouds. 



tion. Indeed slow rotators do not produce GRBs. Very fast 
rotators may not produce GRBs. As can be seen in Table 
the final angular momentum contained in the model 



with Z = Z*gMC! an initial mass of 60 Mq and an initial 
rotational velocity of 500 km s^ 1 is not sufficient. Magnetic 
braking may also reduce the mass range for GRB produc- 
tion. However, since the WO stars spend little or no time 
during the RSG phase and have shorter lifetimes than 
stars around 15 M© stars, the magnetic braking may not 
be as strong for WO as expected by recent calculations 
bv lHeeer et all l|2005h . Finally, these rates have to be con- 
volved in metallicity to compare with observations. The 
T^grb(WO) rates are anyway in much better agreement 
with observation than the rates from all WR stars. The 
very interesting point about WO stars as GRB progen- 
itors is their metallicity dependence. Indeed there is no 
WO predicted at solar or higher metallicities. At very low 
metallicities, the existence of WO stars depends on the 
initial rotational velocity of the stars. If extremely low 
metallicity stars have the same initial surface velocity as 
solar metallicity stars, no WO stars can form at extremely 
low metallicity. However, if the same amount of angular 
momentum as solar metallicity stars is contained in ex- 
tremely low metallicity stars, WO stars may form (see 
Meynet, Ekstrom and Maeder, in prep.). Observations of 
host galaxies of GRB show that these galaxies are sub- 
luminous and blue, pointing in the dire ction that G RB 
produ ction is favoured at low metallicity I Le Floc'h et alJ 
2003J) as is predicted in this study. 

5. Conclusion 

We present the evolution of rotation in models of massive 
rotating single stars covering a wide range in mass and 
metallicity. These models reproduce very well observations 
during the early stages of the evolution (especially WR 
populations and ratio be tween type II and type Ib,c at 
different metallicities, see lMevnet fc Maederll2005l) . 

Our models predict the production of fast rotating 
black holes. Models with large initial masses or high 
metallicity end with less angular momentum in their cen- 
tral remnant with respect to the break-up limit for the 
remnant. Many WR star models satisfy the three main 
criteria (black hole formation, loss of hydrogen-rich en- 
velope and enough angular momentum to form an ac- 
cretion disk around the black hole) for long and soft 
gamma -ray bursts (G RB) production via the collapsar 
model l|Wooslevl |l993) . Considering all types of WR stars 
as GRB progenitors, we predict too many GRBs com- 
pared to observations. If we consider only WO stars (type 
Ic supernovae as is the case for SN2003dh/GRB030329 
iMazzali et"ai]l2003l) as GRB progenitors, the GRB pro- 
duction rates are in much better agreement with obser- 
vations. WO stars are produced only at low metallicities 
in the present grid of models. Although the numbers esti- 
mated in this work are to be taken with caution, the trend 
of the effects found is clear. The prediction for the metal- 
licity dependence can be tested by future observ ati ons. 

Some recent stu d ies llHeeer et alJ 120051: 
IPodsiadlowski et alJ Hil iTzza.rd et alJ l20nl suggest 
that single stars can not retain enough angular mo- 
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mentum for GRBs production and predict that binary 
stars a re the progenitors of GRBs. However. Ilzzard et alJ 
()2004[) use a simple description for angular momentum 
and their model without magnetic b raking are not con - 
sistent with our mod els and those of iHeeer et al.1 l|2000|) . 
iFrver fe Hegerl l|2005j) studied massive binary star mergers 
in detail. They found that, in general, the final angular 
momentum contained in the central remna nt is similar in 
merge d stars and in single massive stars. IPetrovic et alJ 
( 2005) find t hat, inc l uding the effects of magnetic fields 
according to ISpruitl l|2002ft . neither the single nor the 
binary star models that they calculated retain enough 
angular momentum to form GRBs. Finally, magnetic 
braking which is implemented mainly to reproduce the 
observed pulsar periods (which is a very important 
observation), may be much less efficient for WO stars 
than for NS forming stars because WO stars spend little 
or no time in the RSG stage and have shorter lifetimes. 

In conclusion, in the present models, the best candi- 
dates for GRB progenitors are, for single massive stars, 
Wolf-Rayet stars of the type WO, which are found only 
at low metallicities (Z~Zsmc) and have masses, M > 50 
Mq. In order to support this prediction, additional ob- 
servations are needed to confirm the link between type Ic 
supernovae and long-soft GRBs. 
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Table 4. Z=0.004 (SMC): angular and momenta of inertia for rotating stellar models at various stages of their 
evolution. For each model, we give the initial mass and velocity a s well as the remnant mass estimated from the 
carbon-oxygen (CO) core mass using the relation from lMaederl l)l992(l . We used the CO core at the end of the evolution 
for the calculation: at the end of core He-burning in general and at the end of Si-burning for the models from paper 
XII. Then, column 1 is the evolutionary stage to which the values correspond, column 2 is the total mass at the given 
stage, column 3 is angular momentum contained in the remnant, C lcm — Jq 1 '"" 1 2/3 £1 r 2 dm in units of 10 50 g cm 2 , 

column 4 is the moment of inertia of the remnant, T lcm — J^ 1 '"" 1 r 2 dm in units of 10 55 gem 2 , column 5 is the specific 
angular momentum at the remnant edge, j re m, in units of 10 16 cm 2 s , column 6 is the average specific angular 
momentum in the remnant, j rcm = C Icm /M Icm , in units of 10 16 cm 2 s _1 , column 7 is the specific angular momentum 
at the mass coordinate 1.56 Mq, ji.56, in units of 10 16 cm 2 s , column 8 is the angular momentum contained in the 
inner 1.56 M , £1.56 = Jq 6 2/3flr 2 dm in units of 10 50 gem 2 s _1 , column 9 is the specific angular momentum at 
the mass coordinate 2.5 Mq, J2.5, in units of 10 16 cm 2 s _1 , column 10 is the angular momentum contained in the inner 
2.5 M Q , £ 2 .5 = Jq 5 ° 2/357 r 2 dm in units of 10 50 g cm 2 s^ 1 . Assuming that a neutron star with a baryonic mass, 
Mb = 1.56M , and with a radius, R = 12 km, would form from the models, we calculated the two following quantities: 
(column 11) the ratio of the NS angular velocity to its Keplerian angular velocity, H/Hk (NS). (column 12) the initial 
rotation period of the neutron star, V m t in units of milli-seconds. 



stage 


Mass 


£rom 


Item 


jrcm 


Jrcm 


jl.56 


-Cl.56 


j'2.5 


£2.5 


(n/n K ,Ns) 








10 50 


10 55 


10 16 


10 16 


10 16 


10 50 


10 16 


10 50 










[Mq] 


[££ 


^1 


[gem 2 ] 




rf] 




r—"i 




r—"i 






M 










30 M e 


! Vi n i 


= 300 kins 


~\M xem 


= 3.73 M 













ZAMS 


30.0 


11 


.970 


1.8E+00 


27.13 


16.14 


14.24 


2.875 


20.07 


6.190 


17 


20 


0.035 


end H 


28.7 


4. 


309 


4.6E-01 


10.06 


5.81 


5.03 


1.003 


7.24 


2.190 


6. 


00 


0.101 


start He 


28.5 


4 


.123 


9.3E-02 


9.88 


5.56 


4.73 


0.934 


6.94 


2.063 


5. 


59 


0.108 


end He 


18.9 


3 


.153 


2.8E-02 


7.50 


4.25 


3.64 


0.683 


5.30 


1.534 


4. 


on 


0.148 










40 M F 


, Uini 


= 300 kms 


-\M rem 


= 5.40 M 













ZAMS 


40.0 


21. 


.882 


3.8E+00 


34.63 


20.39 


13.86 


2.789 


19.45 


6.164 


16 


6!) 


0.036 


end H 


36.8 


6 


805 


1.3E+00 


11.05 


6.34 


4.20 


0.839 


5.97 


1.870 


5. 


02 


0.121 


start He 


36.5 


6 


367 


2.0E-01 


10.59 


5.93 


3.84 


0.760 


5.52 


1.709 


4. 


.55 


0.133 


start WR 


27.0 


6 


.012 


1.7E-01 


9.97 


5.60 


3.64 


0.721 


5.22 


1.552 


4. 


32 


0.140 


end He 


22.3 


5. 


176 


5.4E-02 


8.58 


4.82 


3.13 


0.592 


4.49 


1.337 


3. 


54 


0.171 










60 M G 


, Uini 


= 300 kms 


-\M ram 


= 8.52 M 













ZAMS 


60.0 


44 


629 


9.7E+00 


44.97 


26.35 


12.95 


2.743 


18.08 


5.775 


16 


41 


0.037 


end H 


51.8 


11 


.797 


3.2E+00 


12.16 


6.97 


3.32 


0.700 


4.67 


1.482 


4. 


19 


0.145 


start He 


50.9 


10 


968 


5.4E-01 


11.52 


6.48 


3.01 


0.632 


4.26 


1.344 


3 


.78 


0.160 


start WR 


44.2 


10 


697 


4.9E-01 


11.23 


6.32 


2.94 


0.617 


4.16 


1.312 


3 


69 


0.164 


end He 


28.5 


8. 


.076 


6.2E-02 


8.70 


4.77 


2.21 


0.431 


3.14 


0.945 


2. 


58 


0.235 










60 M e 


, Uini 


= 500 kms 


"\AW, 


= 3.79 M 













ZAMS 


60.0 


18 


.372 


2.5E+00 


38.32 


24.36 


20.33 


4.306 


28.38 


9.066 


25. 


.76 


0.024 


start WR 


50.3 


4 


.734 


2.1E+00 


9.96 


6.28 


5.21 


1.098 


7.32 


2.322 


6. 


57 


0.092 


end H 


36.7 


1. 


920 


8.0E-01 


4.23 


2.55 


2.19 


0.461 


3.09 


0.978 


2. 


76 


0.219 


start He 


36.3 


1. 


.754 


1.3E-01 


3.90 


2.33 


1.99 


0.418 


2.83 


0.890 


2. 


50 


0.242 


end He 


12.4 


1. 


.102 


2.8E-02 


2.63 


1.46 


1.24 


0.228 


1.82 


0.518 


1. 


37 


0.443 










120 Ifn, Uini 


= 300 km s 


-\M rOT 


, = 5.14 M 










ZAMS 


119.9 


16 


778 


5.4E+00 


25.59 


16.40 


11.09 


2.580 


15.37 


5.408 


15. 


44 


0.039 


start WR 


97.7 


5. 


849 


4.6E+00 


8.96 


5.72 


3.85 


0.893 


5.35 


1.876 


5. 


34 


0.113 


end H 


53.0 


1. 


065 


1.3E+00 


1.73 


1.04 


0.72 


0.154 


1.01 


0.323 


0. 


92 


0.655 


start He 


52.5 


1. 


052 


2.5E-01 


1.71 


1.03 


0.71 


0.153 


1.00 


0.319 


0. 


91 


0.664 


end He 


17.2 


0. 


478 


5.3E-02 


0.84 


0.47 


0.31 


0.059 


0.45 


0.131 


0. 


35 


1.724 
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Table 5. Z=0.008 (LMC): angular and momenta of inertia for rotating stellar models at various stages of their 
evolution (sec Table ^ for the legend). 



stage 


Mass 


£ [cm 


Txem 


jrem 


J'rcm 


il.66 


•Cl.56 


32.5 


£2.5 


(n/n K , NS ) 


("ProO 






10 50 


10 55 


10 16 


10 16 


10 16 


1Q 50 


10 16 


10 50 








f Mr,] 
L JVJ W J 


I s J 


r 2i 

[g cm j 


L s J 


L a J 


l s J 


L s J 


.cm*, 

L s J 


[««»», 

L s J 




r I 

[msj 












— OUU Kill s 


- 1 M 


2 73 ]VI 


3 








ZAMS 


30.0 


12.106 


1.9E+00 


27.46 


16.34 


14.42 


2.911 


20.33 


6.268 


17.42 


0.035 


end H 


27 5 


4 072 


5 3E-01 


9 51 


5 49 


4.76 


949 


6 85 


2 073 


5 68 


107 


start tic 


27.3 


3.880 


9.7E-02 


9.30 


5.24 


4.46 


0.880 


6.54 


1.944 


5.27 


0.115 


start WR 


16.3 


3.365 


7.9E-02 


8.13 


4.54 


3.96 


0.744 


5.77 


1.670 


4.45 


0. 136 


end tic 


12.0 


2.804 


1. 1E-02 


6.86 


3.78 


3.27 


0.597 


4.81 


1.381 


3.57 


0. 170 








4U 1V1 (7) 




— OUU K.II1 o 


-i M rBm 


— 5 21 M 










ZAMS 


40.0 


20.719 


3.9E+00 


33.88 


19.99 


13.93 


2.803 


19.54 


6.193 


16.77 


0.036 


end H 


35.3 


5.935 


1.2E+00 


9.95 


5.73 


3.90 


0.777 


5.54 


1.733 


4.65 


0. 130 


start He 


34.8 


5.510 


2.0E-01 


9.45 


5.32 


3.54 


0.701 


5.08 


1.575 


4. 19 


0. 144 


start WR 


28.1 


5.338 


1.7E-01 


9.14 


5.15 


3.43 


0.680 


4.93 


1.464 


4.07 


0. 149 


end He 


17.2 


4.052 


6.5E-03 


7.11 


3.91 


2.62 


0.493 


3.78 


1.118 


2.95 


0.205 








60 M 


i ^ini 


= 300 kms 




= 4.91 M 











ZAMS 


60.0 


17.965 


4.1E+00 


29.67 


18.40 


13.02 


2.757 


18.17 


5.804 


16.49 


0.037 


end H 


48.4 


5.156 


7.7E-01 


8.64 


5.28 


3.68 


0.775 


5.18 


1.642 


4.64 


0.131 


start He 


46.9 


4.594 


2.1E-01 


8.02 


4.71 


3.36 


0.705 


4.75 


1.499 


4.22 


0.143 


start WR 


43.8 


4.523 


1.8E-01 


7.90 


4.63 


3.31 


0.694 


4.68 


1.476 


4.16 


0.146 


end He 


16.4 


2.688 


4.1E-02 


4.98 


2.75 


1.92 


0.358 


2.78 


0.804 


2.14 


0.283 








120 M B , v ini 


= 300 km s 


~\A/ rc „ 


, = 4.06 M 








ZAMS 


119.9 


11.385 


4.0E+00 


21.45 


14.11 


11.00 


2.560 


15.25 


5.366 


15.32 


0.040 


start WR 


93.8 


3.394 


3.5E+00 


6.42 


4.21 


3.27 


0.759 


4.55 


1.595 


4.54 


0.133 


end H 


36.6 


0.525 


7.5E-01 


1.07 


0.65 


0.53 


0.106 


0.74 


0.235 


0.63 


0.958 


start He 


36.3 


0.519 


1.6E-01 


1.06 


0.64 


0.52 


0.104 


0.73 


0.232 


0.62 


0.970 


end He 


13.3 


0.277 


6.3E-03 


0.62 


0.34 


0.28 


0.052 


0.41 


0.116 


0.31 


1.950 
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Table 6. Z=0.020 (solar): angular and momenta of inertia for rotating stellar models at various stages of their evolution 
(see Table 0] for the legend). 



stage 


Mass 


£,„ 


x rcm 


jrem 


3 rem 


il.66 


•Cl.56 


32.5 


£2.5 


(n/n K] Ns) 


CPrat) 






10 50 


10 55 


10 16 


10 16 


10 16 


10 50 


10 16 


10 50 








f Mr,] 
l lyA J 


[ — 


— ] 


r 2i 

IS cm j 


L s J 


r cm2 i 


I s J 


f- "1 




L s J 




r -I 
[I11SJ 










Q A I ^ 

a mq. 




— OUU Kill 3 




= 1.18 Mq 








ZAMS 


9.0 


1 


.632 


1.5E-01 


11.51 


6.95 


14.30 


2.592 


21.10 


5.924 


15.51 


0.039 


end H 


8 8 





.744 


7 1E-02 


5 49 


3 17 


7.07 


1.196 


11.90 


2 976 


7 16 


085 


start tic 


8.8 





.719 


4.5E-03 


5.38 


3.06 


7.47 


1.189 


0.71 


2.806 


7. 12 


0.085 


end tic 


8.4 





.639 


2.1E-03 


5.14 


2.72 


5.75 


1.057 


13.14 


2.640 


6.32 


0.096 










12 Mq 




— ^nn L-m c" 

— OUU Kill 3 


-i m 


= 1.46 M 










ZAMS 


12.0 


2 


.446 


2.7E-01 


13.86 


8.42 


14.57 


2.730 


21.15 


6.063 


16.33 


0.037 


end H 


11.5 


1 


.025 


1.1E-01 


6.05 


3.53 


6.40 


1.149 


9.92 


2.642 


6.88 


0.088 


start He 


11.5 


1 


.021 


9.4E-03 


5.98 


3.51 


6.33 


1.144 


10.23 


2.635 


6.85 


0.088 


end He 


10.3 





.863 


3 8E-03 


5.46 


2.97 


5.84 


0.975 


8.33 


2.325 


5.83 


0. 104 


start O 


10.2 





.688 


4.2E-05 


5.53 


2.37 


6.86 


0.821 


8.63 


2.141 


4.91 


0. 123 










15 A4q 




— ^nn L-m c" 

— OUU Kill 3 


-i m 


= 1.85 M 


© 








ZAMS 


15.0 


3 


.682 


4.5E-01 


16.79 


10.01 


14.76 


2.762 


21.21 


6.196 


16.52 


0.037 


end H 


14.2 


1 


.455 


2.0E-01 


6.87 


3.95 


5.95 


1.080 


9.02 


2.506 


6.46 


0.094 


start He 


14.2 


1 


.452 


1.8E-02 


6.87 


3.95 


6.19 


1.073 


8.96 


2.500 


6.42 


0.094 


end He 


10.4 


1 


.182 




5.78 


3.21 


4.93 


0.872 


7.88 


2.065 


5.22 


0. 116 




10.3 





.041 


1.9E-05 


4.47 


2.56 


4.10 


0.729 


11.40 


2.015 


4.36 


0. 139 










£\J IVl 




— OUU Kill S 


-1 ]y£ 


= 2.57 M 










ZAMS 


20.0 


6 


.590 


9.3E-01 


21.73 


12.91 


14.95 


2.896 


°21.30 


6.304 


17.33 


0.035 


end H 


18.2 


2 


.322 


3.5E-01 


7.92 


4.55 


5.25 


0.997 


7.74 


2.217 


5.96 


0.102 


start He 


17.5 


2 


.208 


3.7E-02 


7.81 


4.33 


4.97 


0.924 


7.61 


2.105 


5.53 


0. 109 


end He 


8.8 


1 


.781 


1.3E-02 


6.31 


3.49 


4.08 


0.740 


6.15 


1.698 


4.43 


0. 137 




8.8 


1 


.652 


6 6E-05 


6.61 


3.24 


2.80 


0.627 


6.50 


1.565 


3.75 


0. 161 










ZD JW0 




— OUU Kill S 


-i M 


= 3.06 M 











Z A 


25 


8 


.990 




24 85 


14 79 


15 06 


2.976 


21.33 


6 403 


17 81 


034 


end H 


21.8 


2 


.898 


5.1E-01 


8.24 


4.77 


4.80 


0.935 


6.97 


2.045 


5.60 


0. 108 


start He 


21.6 


2 


.786 


6.5E-02 


8.19 


4.58 


4.55 


0.873 


6.80 


1.944 


5.22 


0.116 


start WR 


13.5 


2 


.457 


5.3E-02 


7.26 


4.04 


4.11 


0.759 


6.07 


1.714 


4.54 


0. 133 


end He 


10.2 


2 


.143 


2.1E-02 


6.38 


3.53 


3.57 


0.658 


5.31 


1.492 


3.94 


0. 154 


end Si 


10.0 


1 


.644 


7.1E-05 


3.72 


2.71 


1.90 


0.600 


3.20 


1.186 


3.59 


0. 169 










An A/7" 




QflD L-m c ' 

— OUU Kill S 


-i m 


= 3.85 M 











Z A 

Zj t\. 1V± O 


40 


12 


.549 




26 82 


16 38 


13 88 


2.793 


19.48 


6 173 


16 71 


036 


end H 


32.9 


3 


.035 


8.9E-01 


6.61 


3.96 


3.32 


0.662 


4.71 


1.476 


3.96 


0.153 


start He 


32 3 


2 


.847 


1 3E-01 


6 30 


3 72 


3 07 


0.609 


4.42 


1 348 


3 65 


166 




23 9 


2 


.659 


1 0E-01 


6 02 


3.47 


2 95 


0.586 


4.24 


1 260 


3 51 


173 


end He 


12.7 


1 


.976 


3.0E-02 


4.68 


2.58 


2.18 


0.407 


3.19 


0.910 


2.44 


0.249 




12.6 


1 


.399 


2.9E-04 


3.00 


1.83 


1.86 


0.361 


0.98 


0.849 


2. 16 


0.281 










60 M e 


i l-'ini 


— ^nn L-m 

— OUU ivlll o 


- 1 m 


= 4.32 M 











ZAMS 


60.0 


14 


.451 


3.8E+00 


26.77 


16.81 


12.90 


2.732 


18.01 


5.753 


16.35 


0.037 


start WR 


44.4 


3 


.644 


3.1E+00 


7.08 


4.24 


3.34 


0.705 


4.70 


1.492 


4.22 


0. 144 


end H 


31.4 


2 


.548 


7.2E-01 


5.01 


2.96 


2.32 


0.454 


3.28 


0.991 


2.71 


0.223 


end He 


14.6 


1 


.662 


2.0E-02 


3.49 


1.93 


1.49 


0.276 


2.16 


0.622 


1.65 


0.367 


end Si 


14.6 


1 


.193 


5.0E-04 


1.71 


1.39 


1.00 


0.250 


2.17 


0.615 


1.50 


0.404 










85 Mq 


l ^ini 


= 300 kins 


~\M xem 


= 3.78 M 











ZAMS 


84.9 


10 


.546 


3.4E+00 


21.97 


14.04 


11.79 


2.638 


16.39 


5.461 


15.78 


0.038 


start WR 


64.5 


4 


.401 


3.1E+00 


9.22 


5.86 


4.90 


1.094 


6.84 


2.271 


6.54 


0.093 


end H 


23.6 





.727 


5.6E-01 


1.67 


0.97 


0.85 


0.161 


1.21 


0.357 


0.97 


0.627 


end He 


12.3 





.528 


1.4E-02 


1.27 


0.70 


0.60 


0.110 


0.88 


0.253 


0.66 


0.917 










120 Mq, u ini 


= 300 km s 


" 1 ,M rci; 


, = 3.54 M B 








ZAMS 


119.9 


8 


.802 


3.5E+00 


18.68 


12.50 


10.55 


2.455 


14.63 


5.146 


14.69 


0.041 


start WR 


87.2 


2 


.344 


3.2E+00 


4.99 


3.33 


2.80 


0.650 


3.90 


1.367 


3.89 


0.156 


end H 


26.2 





.386 


6.1E-01 


0.93 


0.55 


0.50 


0.099 


0.71 


0.214 


0.59 


1.026 


start He 


26.0 





.376 


1.2E-01 


0.91 


0.53 


0.49 


0.096 


0.69 


0.208 


0.57 


1.055 


end He 


11.3 





.248 


1.6E-02 


0.64 


0.35 


0.32 


0.059 


0.47 


0.133 


0.35 


1.719 
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R. Hirschi et al.: Stellar evolution with rotation XIII: 



Table 7. Z=0.040 (GC): angular and momenta of inertia for rotating stellar models at various stages of their evolution 
(see Table 0] for the legend). 



stage 


Mass 


L 




x rcm 


jrem 


Jrcm 




J1.66 


•Cl.56 




32,5 


£2.5 


(n/n K , NS ) 


CProt) 






1Q 50 


10 55 


10 16 


10 16 




10 16 


10 50 




10 16 


10 50 












f Mr,] 
L JVJ 3 J 


[££ 


— ] 


r_.__.2i 

[g cm j 


L s J 


n 2 ] 




I s J 


[«<»», 
L s -1 




r cm2 i 


[ — 


— 1 








ms| 










20 Mq 




-- 300 kins" 


- 1 M 


— 2 68 


M 


3 














ZAMS 


20.0 


7 


.033 


1.1E+00 


22.27 


13.19 




14.82 


2.869 




21.11 


6. 


245 


17 


.17 





.035 


end H 


16 9 


2 


.223 


3 0E-01 


7 29 


4.17 




4 66 


882 




6 87 


1. 


.966 


5 


.28 





.115 


start tic 


16.8 


2 


.176 


4.8E-02 


7.43 


4.08 




4.52 


0.840 




6.93 


1. 


.915 


5 


.03 





.120 


end tic 


9.2 


1 


.742 


1.1E-02 


5.92 


3.27 




3.69 


0.670 




5.55 


1. 


.535 


4 


.01 





.151 










25 Mq 




= 300 kins" 


_1 M 


= 3 13 


M 


3 














ZAMS 


25.0 


9 


.300 


1.6E+00 


25.26 


14.96 




15.04 


2.977 




21.31 


6. 


396 


17 


.81 





.034 


end H 


21.3 


1 


.980 


3.2E-01 


5.48 


3.18 




3.18 


0.623 




4.58 


1. 


.352 


3 


.73 





.163 


start He 


20.9 


1 


.872 


8.4E-02 


5.24 


3.01 




2.99 


0.583 




4.35 


1. 


.273 


3 


.49 





.174 


start WR 


18.5 


1 


.854 


6.6E-02 


5.18 


2.98 




2.96 


0.554 




4.30 


1. 


.249 


3 


.31 





.183 


end He 


9.6 


1 


.471 


1.9E-02 


4.31 


2.37 




2.35 


0.424 




3.51 





.976 


2 


.54 





.239 










40 Mq- 




= 300 kins 


- 1 M reix 


— 3,57 


M 
















ZAMS 


40.0 


10 


.867 


2.5E+00 


24.86 


15.30 




13.61 


2.737 




19.10 


6 


051 


10 


.38 





.037 


start WR 


31.1 


2 


.400 


1.7E+00 


5.74 


3.38 




3.07 


0.612 




4.35 


1. 


.308 


3 


.66 





.165 


end H 


29.9 


2 


.343 


5.5E-01 


5.62 


3.30 




2.99 


0.596 




4.25 


1. 


.275 


3 


.57 





.170 


start He 


29.7 


2 


.170 


1.2E-01 


5.27 


3.05 




2.75 


0.545 




3.95 


1. 


.172 


3 


.26 





.186 


end He 


11.4 


1 


.516 


2.3E-02 


3.87 


2.13 




1.91 


0.356 




2.81 





.798 


2 


.13 





.285 












— 300 


kms -1 


,M, 


■cm — 3 


02 Mq 


M 


WR ( Z) 














ZAMS 


40.0 


8 


.250 


i QF,-i_nn 
± . ,j ____ \j \ i 


21.95 


13.72 




13.61 


2.737 




19.10 


6 


051 


16 


.38 





.037 


start WR 


31.1 


1 


.804 


i 3E-i-on 


5.03 


3.00 




3.07 


0.612 




4.35 


1. 


308 


3 


.66 





.165 


end H 


29.5 


1 


.754 


4.0E-01 


4.91 


2.92 




2.98 


0.594 




4.24 


1. 


.271 


3 


.55 





.170 


start He 


29.2 


1 


.623 


8.9E-02 


4.59 


2.70 




2.75 


0.544 




3.94 


1. 


.172 


3 


.26 





.186 


end He 


9.0 


1 


.087 


1.7E-02 


3.29 


1.81 




1.84 


0.331 




2.76 





.770 


1 


.98 





.305 










fid Kfr> D- ■ 


— 300 


km s 1 


,M S 


■cm — 1 


.94 Mq , 


M 


( 7\ 
WR \ ^ J 














start WR 


44.9 





.971 


9.9E-01 


3.83 


2.52 




3.27 


0.691 




4.60 


1. 


461 


4 


.13 





.146 


end H 


9.3 





.169 


1.3E-01 


0.76 


0.44 




0.63 


0.116 




0.94 





.265 





.69 





.872 


start He 


9.1 





.166 


3.0E-02 


0.75 


0.43 




0.62 


0.114 




0.93 





.261 





.68 





.889 


end He 


4.8 





.136 




0.66 


0.35 




0.53 


0.091 




0.87 





.221 





.54 


1 


.111 














= 300 kins" 


- 1 m 


— 2 57 


M 


3 














ZAMS 


84.8 


5 


.316 


1.9E+00 


15.57 


10.40 




10.98 


2.456 




15.27 


5. 


.086 


14 


.70 





.041 


start WR 


62.6 


2 


.056 


1.8E+00 


6.05 


4.02 




4.24 


0.947 




5.93 


1. 


907 


5 


.67 





.107 


end H 


15.6 





.334 


2.5E-01 


1.12 


0.65 




0.75 


0.142 




1.09 





.318 





.85 





.711 


start He 


15.3 





.329 


5.7E-02 


1.11 


0.64 




0.74 


0.140 




1.09 





314 





.84 





.725 


end He 


7.3 





.237 


1.2E-02 


0.85 


0.46 




0.54 


0.096 




0.83 





.225 





.58 


1 


.051 










85 M , v iDi 


= 300 


kms" 1 


,M, 


— 1 


95 Mq 


M 


wrA^/ 














ZAMS 


84.8 


3 


.455 


1.3E+00 


12.84 


8.90 




10.98 


2.456 




15.27 


5. 


.086 


14 


.70 





.041 


start WR 


62.6 


1 


.334 


1.2E+00 


4.97 


3.44 




4.24 


0.947 




5.93 


1. 


.967 


5 


.67 





.107 


end H 


11.6 





.202 


1.5E-01 


0.90 


0.52 




0.75 


0.137 




1.11 





.313 





.82 





.737 


start He 


11.3 





.204 


3.0E-02 


0.92 


0.53 




0.76 


0.138 




1.14 





.317 





.82 





.735 


end He 


4.7 





.155 


7.2E-03 


0.75 


0.40 




0.60 


0.102 




0.98 





.249 





.61 





.993 










120 M ra , v ini ■- 


= 300 km s 


~\M rcr 


n = 2.53 


M 


-3 














ZAMS 


118.8 


5 


.616 


2.2E+00 


15.80 


11.15 




11.29 


2.626 




15.65 


5. 


.505 


15 


.71 





.039 


start WR 


86.5 


1 


.314 


2.1E+00 


3.71 


2.61 




2.64 


0.613 




3.67 


1. 


.287 


3 


.06 





.165 


end H 


14.0 





.132 


3.1E-01 


0.45 


0.26 




0.31 


0.058 




0.45 





.129 





.35 


1 


.755 


start He 


13.8 





.130 


5.7E-02 


0.45 


0.26 




0.30 


0.057 




0.44 





.127 





.34 


1 


.784 


end He 


7.1 





.103 


1.3E-02 


0.38 


0.21 




0.24 


0.043 




0.37 





101 





.26 


2 


.341 










120 M Q , v in 


i = 300 


1 kms -1 


,M rem = 1 


.96 M 0! 


M-wr(Z) 














ZAMS 


119.6 


3 


.969 


1.5E+00 


13.68 


10.16 




11.68 


2.717 




16.20 


5. 


.697 


16 


.26 





.037 


start WR 


86.5 





.896 


1.4E+00 


3.10 


2.29 




2.64 


0.613 




3.67 


1. 


.288 


3 


.67 





.165 


end H 


9.6 





.083 


1.6E-01 


0.37 


0.21 




0.30 


0.055 




0.45 





.127 





.33 


1. 


.842 


start He 


9.4 





.082 


3.1E-02 


0.36 


0.21 




0.30 


0.054 




0.45 





.126 





.32 


1 


.875 


end He 


4.8 





.068 


6.8E-03 


0.32 


0.17 




0.26 


0.044 




0.42 





.108 





.26 


2 


.296 



